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Abstract. We find a general solution to the unique 7th order ODE admitting ten dimen- 
sional group of contact symmetries. The integral curves of this ODE are rational contact 
curves in P 3 which give rise to rational plane curves of degree six. The moduli space of 
these curves is a real form of the homogeneous space Sp(4)/SL(2). 



The 7th order ODE 

10(y( 3 ))y 7 )-70(y( 3 ))VV^ = 0, (1) 

where y = y(x) and y( k ' = d k y/dx k has recently appeared explicitly [121 El HI E] or 
implicitly [I] in several different contexts. It is the unique (up to contact transformations) 
equation admitting ten-dimensional algebra of contact symmetries [12] , and the aim of this 
note is to show that its general solution is given by a degree six rational curve of the form 

y 3 + a(x)y 2 + (3(x)y + j(x) = 0, (2) 

where (a,/3, 7) are a quadratic, a quartic, and a sextic respectively with the coefficients 
depending on seven parameters as in formula f)13|) . 

In fact the symmetry algebra of ([I]) was known to Lie [9] who also proved that this is the 
maximal algebra of contact vector fields on the plane. It is quite possible that equation (JTJ) 
and its general solution were also known to Lie. We have been unable to find it in any of 
Lie's works , but we would be grateful to hear from anyone who has (the earliest reference 
to the equation - but not its solution - we are aware of is [10J). 

Contact Lie algebras. Let U C I 2 be an open set and let F(T*U) be a projectivised 
cotangent bundle with a contact one-form u. A curve 7 C ¥(TU) is called contact if 
u)\y = 0. A contact transformation is a map / : ¥(TU) — > ¥(TU) which takes contact 
curves into contact curves. Equivalently f*(ui) = \ui for some function A. Let (x,y) be 
the local coordinates on U and let z parametrise the fibers of ¥(TU), so that we can set 
u) = dy — zdx. Consider a one parameter group of contact transformations. Close to the 
identity, this is characterised by a contact vector field X such that the contact condition 
Cx^ = cw holds, where Cx = d(X J ) + X J d is the Lie derivative and c is some function. 
The contact condition implies that locally there exist a function H = H(x, y, z) such that 

x H = -(d z H)d x + {h- zd z H)d y + (d x H + zd y H)d z . (3) 

If H = a(x, y) + zb(x, y) then Xh generates a prolongation of a family of point transforma- 
tions f : U — > U. Otherwise it generates a proper contact flow. 

The symmetry. A remarkable result of Lie is that a maximum dimension of a Lie algebra 
of proper contact vector fields on the plane is ten. This maximal, ten-dimensional Lie 
algebra is generated by vector fields Q corresponding to functions 

1, x, x 2 ,y, z, xz, x 2 z — 2xy, z 2 , 2yz — xz 2 ,4xyz — Ay 2 — x 2 z 2 . (4) 

This algebra is isomorphic to so (5), or equivalently to sp(4) (this is, up to a choice of the 
real form, the maximal subalgebra of the 11 dimensional symmetry algebra of the trivial 
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7th order ODE). The first seven generators correspond to prolonged point vector fields and 
the last three are the proper contact vector fields. 

The algebra generates all contact symmetries of equation ([T]): if y(x) is a solution to 
([1]) then so is y(x), where 

~ -/ \ ~ -/ x ~ ~, x d x y + zd y y 

x = x{x,y,z,Ci), y = y{x,y,z,Ci), z = z(x,y, z,Ci) = — — — d z y = zd z x 

O x X + ZOyX 

is the contact transformation generated by the vector fields and c%, . . . , cio are parame- 
ters of this transformation. Equivalently, any of the generators (dJ satisfies the linearisation 
of ([T]) when z = y' . 

Another result of Lie [S] is that a maximal dimension of the contact symmetry algebra of 
an ODE of order n > 3 is (n + 4) , with maximal symmetry occurring if only if the ODE is 
contact equivalent to a trivial equation ?/ n ) = 0. Therefore equation ([1]) is of submaximal 
type [IT] - it is not equivalent to the trivial equation and its symmetry algebra has the 
largest possible dimension. Up to the contact equivalence ([1]) is the unique 7th order ODE 
with this property. 

The solution. We verify that the algebraic curve y 2 + x{x — l) 3 = solves ([I]). This curve 
is parametrised by 

x{t) = wrv v{t) = -wrw (5) 

The 7-dimensional subalgebra of dH) consisting of prolonged point symmetries integrates 
to 

y -> c A y + ci + c 2 x + c 3 x 2 , x -> c 5 x + c 6 , and y -> — — ^, x -> — ^ — • (6) 

(1 + C-jxY 1 + CjX 

The contact transformations generated by ^ with H$ = z 2 ,Hq = 2yz — xz 2 and H\q = 
Axyz — \y 2 — x 2 z 2 respectively are given by 

x = x - 2c s z, y = y - c 8 z 2 , z = z (7) 

x(l + c 9 z) - 2c 9 y . y(l - 2c 9 z) + c 9 xz 2 z 

x = : , y= 7z To , z=- , and 

1 — c 9 z (1 — c 9 z) z 1 — C 9 Z 

x y + 4ci ?/ 2 - ^c w xyz + c 1Q x 2 z 2 



1 + 4ci y - 2c w xz ' (1 + Ac 10 y - 2c 10 xz) 2 1 + Ac w y - 2c 10 xz 

where eg, eg, ciq G R. 

Applying the group of point transformations to the given solution yields the six parameter 
family of solutions to dU) given by a family of algebraic curves of degree four 

(y + Q) 2 + P = (8) 

where Q = Q{x) is an arbitrary quadratic, and P = P(x) is a quartic with one simple and 
one triple root. 

Recall [13] that a point p on an algebraic curve f(x,y) = in CP 2 is singular if the 
partial derivatives f x and f y vanish at p. Moreover p has multiplicity m if all (m — l)st 
derivatives of / vanish at p but at least one mth derivative does not. A singular point is 
called ordinary if the tangents to all branches at the point are distinct. Any singular point 
characterised by a triple of integers (m,5,r), where m is the multiplicity, r is the number 
of branches and 5 is the number of multiplicity two ordinary singular points concentrating 
at p. The arithmetic genus of a curve is given by 



g 
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where d is the degree of the curve, and the summation is taken over all singular points. 
The curves in the family ([5]) have two singularities: a cusp of type (2,1,1) at (x,y) = 
(xq, —Q(xo)), where xo is the triple root of P and a point at oo of type (2, 2, 2). Calculating 
the genus of curves in this family yields 

& 2 

Therefore the family is rational. Applying the point transformations in ([6]) to the parametri- 
sation ([5]) and redefining the constants c\, . . . , cq we find that the rational parametrisation 
is given by 

_ h + ht 2 _ b 4 t 4 + b 3 t 2 + bjt + h 

X[t) ~ b + t 2 ' y[t) ~ (b + t 2 ) 2 

The six parameters in ([8]) are algebraic expressions in the seven parameters (po, . . . , 6g) ° ne 
of which is irrelevant and arises only in the parametrisation. 

The family of curves ([8]) can also be obtained applying Hg and the point transformations 
to the trivial solution y = x 2 . It is not the general solution to (pQ) as it depends on six 
parameters rather than seven. To introduce the additional parameter, and construct the 
general solution we use the contact transformation ([7]) generated by ([3]) with H% = z 2 

x(t) = x(t)-2bz(t), y(t) = y{t) - bz{t) 2 , z(t)=z(t), (9) 
y(t) (4&4&0 " 2h)t 3 - 3b 2 t 2 + (2b 3 b - 4&i)t + b 2 b 



where 



z(t) 



x{t)~ 2(b 5 -b b 6 )(b t + t 3 ) 



and (b,bo, . . . ,be) are constant parameters. The relation ([9]) gives a seven-dimensional 
family of rational contact curves in P 3 . The connected component of the symplectic group 
Sp(A) acts on P 3 and preserves the family ([9]). The symmetry group of any fixed rational 
curve in this family is SL(2), and so the seven-dimensional space of solutions to ([1]) is 
the symmetric space Sp{4)/SL(2). In the holomorphic category, a rational curve in P 3 
can be characterised by a normal bundle, which in our case is N = 0(5) © 0(5). The 
contact modification of the Kodaira theorem described in [2] can be applied to deduce that 
we have constructed a complete analytic family of contact curves: an infinitesimal contact 
deformation of any fixed curve in the family ([9|) also belongs to this family. 

Alternatively, the general solution to ([T]) can be given by an implicit relation 

u(x,y,z)=0, v{x,y,z)=0 (10) 

where (x,y,z) are coordinates on an open set in P 3 . Using ([5D and ^ we find 

u = (y + bz 2 + Q(X)) 2 + P(X), v = AP{z + Q'f + {P'f 

where Q = Q(X) is a quadratic, P = P{X) is a quartic with one simple and one triple 
root, and X = x + 2bz. To find an explicit formula for y(x) pick a real root z of the cubic 
v = 0, substitute this in u = and solve the resulting quadratic for y. Alternatively we 
use the resultant to produce a planar curve birationaly equivalent to (jlOp by eliminating z 
between u and v. Recall [7] that a resultant of two polynomials 

u(z) = u + u\z H h u m z m , v(z) = vq + viz H h v n z n 
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Figure 1. Curve (HI]) with 6=1/2 and curve ([8]) with Q = 0, P = x(x — 1) 




is the determinant of the matrix 



f u m U m —l 

Urn ^m—l 



U 



*m— 1 



Uq 



'Jn-l 



V n -1 



\ 



V n V n -i 



The resultant vanishes if u, v admit a common root. To obtain a manageable formula 
consider (|10p with Q = and P = X(X — l) 3 , and find the resultant of u and v. This 
resultant factorises into two terms, each giving a rational curve of degree six. We choose one 
of these two curves (The second curve is not a solution to the ODE ([1]). It arises because 
the expressions in (|10p are squares of actual solutions and contain the term z 2 . Only one 



of the roots satisfies z = y'). 

(646 + 1024 b 3 ) 



((768 6 2 + 16) x 2 - 768 xb 2 + 28 8 6 2 ) y 2 



+ (264 x 2 b - 108 6 3 + 192 x 4 6 - 72 x6 - 384 x 3 b) y 



(11) 



+ (48: 



27 6 2 + 54x6 2 - 16 x s 



27x 2 b 2 



48 x 5 + 16 x 6 ) = 0. 



This curve does not belong to the class (JS]): it has two (2, 1, 1) cusps and one non-ordinary 
(2,2,2) singularity at oo whereas ([H]) has one cusp (apart from the singularity at oo). The 
curve (|lip possesses the following property: its discriminant is a cube of a quartic Q(x) 
in x with two real roots. These roots correspond to the positions of the finite cusps, as in 
Figure [TJ For any such curve ([2]) the cross-ratio p of the roots of Q is an invariant with 
respect to transformations (0). It is possible to verify that for (TTT]) we have p = e, where 
e = exp (^j)- To bring the curve (jlip to a canonical form we first set the coefficient of y 2 
to zero by a transformation of the form y — > y + c\ + Cix + C3X 
Mobius transformations can be used to set Q(x) = x(x — l)(x 
y, yields 

1 



2 After that the remaining 
p). This, up to scaling of 



y 



3y[x 



1 



1 + e 



1 + e 



•!+£) + 



2(x - (1-ieU (x - (l+iej\ (x-(i-iej\ (x-(-i + ie)^ (x - (i + e)) (x-(-i + e)\ 
The discriminant Q(x) of this curve is proportional to x 3 (x — l) 3 (x — e) 3 . 



0. 
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Using this complex canonical form, we easily can find a real one. To do that we find 
four roots of the sextic in x in the complex canonical form such that their cross-ratio is 
— 1. Using the Mobius transformations, we move these roots to ±1 ± i. The result (up to 
scaling of x and y) is given by 

y 3 + 3(3x 4 - 6x 2 - l)y + 12x(3x 4 + 1) = (12) 

which can be parametrised by 

( ) 3(t 2 + 1) ' V[ > 3(t 2 + l) 2 ' 

To get the general solution for the 7th order ODE (pQ) we apply the point transformations 
to (fT2j) . The resulting degree six rational curve of the form ([2]) is given by 

(c 4 y + ci + c 2 x + c 3 x 2 ) 3 (13) 

+3 ^3 (c 5 x + c 6 ) 4 - 6 (c 5 x + c 6 ) 2 (1 - c-jxf - (1 - c 7 x) 4 ^ (c 4 y + c\ + c 2 x + c 3 x 2 ) 

+12 (c 5 x + c 6 ) (^3 (c 5 x + c 6 ) 4 (1 - c 7 x) + (1 - c 7 x) 5 ) = 0. 
The formula ([2]) is the general solution as the generic initial data 

can be chosen arbitrarily by choosing the coefficients c%, . . . , C7 in the solutions. There exist 
additional singular solutions corresponding to submanifolds of the initial data manifold 
which can not be obtained by any choice of coefficients in (|13p . An example of such 
singular solution with a co-dimension one initial data set is the curve ([8]). The leading 
term in this curve is y 2 and this can not arise in (|13p . This type of behaviour should not 
be confused with the singular orbits of the symmetry group. The connected component of 
the group Sp(A) is transitive on the solution space, and thus (fT3"j) and ([H]) are related by a 
contact transformation. This is indeed how we found (|13|) . 

Conclusions. We have found the general solution of the 7th order ODE (pQ). This equation 
is submaximal as its symmetry algebra is ten-dimensional, wheres the algebra of the trivial 
7th order ODE is eleven-dimensional. The general solution to ([I]) is given by the seven 
dimensional orbit of the point transformations ([6]) acting on the canonical solution ()12p . 
There also exists a six-dimensional submanifold in the space of solutions given by degree 
four curves ([8]) . For each contact transformations generated by H$ , Hg and -f/io in there 
exists a point transformation such that the composition of the two fixes ()12p . This gives a 
stabiliser SL(2) of (fT2|) and finally the solution space Sp(4)/SL(2). 

The analogous, submaximal 5th order ODE characterises conies in MP 2 . In the inhomo- 
geneous coordinates (x, y) the five parameter family of conies is 

y 2 = ci x 2 + c 2 xy + c 3 y + c 4 x + c 5 . 

Eliminating the parameters (ci, . . . ,c$) between this equation and its fourth derivatives 
and substituting in the fifth derivative yields the ODE 

9 (y(2) ) 2 y (5) _ 45y (2) y (3) y (4) + 40(y (3) )3 = q 

This construction goes back to Halphen [8] who wrote the equation as ((y( 2 )) _2 / 3 )( 3 ) = 0. 

The seven-dimensional space of solutions M to (HJ carries a GL(2,M) structure in the 
sense of OS]: TM has a pointwise identification with a vector space of homogeneous degree 
six polynomials in two variables. The five Wiinschmann-Doubrov-Wilczynski invariants 
vanish on the ODE (PQ) which implies [3] that the linearisation of ([1]) is equivalent to a 
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trivial ODE 5y^ = 0. Moreover [6], M also admits a conformal structure of signature 
(3,4). The null vectors of this structure correspond to the six order polynomials 

aix 6 + 6a2X 5 + 15a3X 4 + 20a4X 3 + 15asx 2 + Ga^x + 07 

with vanishing quadratic invariant [7j 

a\(i7 — 60206 + 1^0305 — lOa^. 

In [6] it was shown that the conformal structure associated to ([1]) contains a metric with 
weak holonomy G2: there exists a three form on M such that 

d<j> = A * cj), d* 4> = 

where A is a constant, and * is the Hodge operator. Taking an analytic continuation 
of this structure to the Riemannian signature yields the homology seven-sphere M = 
SO(5)/SO(S) with its canonical weak G2 structure originally constructed by Bryant [I]. 
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